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The aim of the paper is to make a comparison between the unitarity condition 
method and the standard version of the unitarity triangle approach by using as pa- 
rameters four independent moduli \Uij\. This choice is motivated by the measurability 
property and leads to a simple criterion for the separation of unistochastic matrices from 
the double stochastic ones, whose fulfillment is the key point of any global fit for the 
CKM entries. In our formulation both the methods are exact, do not depend on any as- 
sumption as the smallness of some parameter and both can be used to global fits in the 
quark and fermion sectors. Monte Carlo simulations show that the separation criterion 
puts very strong conditions requiring a fine tuning of all the CKM matrix elements. 
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1. Introduction 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix is a lively subject in particle 
physics and the determination of its entries that govern all flavor changing transi- 
tions of quarks and leptons in the Standard Model is an important task for both 
experimenters and theorists. The large interest in the subject is also reflected by the 
workshops organized in the last years whose main subject was the CKM matrix 1-4 , 
one of the main goals of these gatherings being global fits to the CKM matrix 
parameters. The usual method for doing that is the so called standard unitarity 
triangle approach that make use of the orthogonality property between the first 
column and the third one of the recommended form for the CKM matrix, see Refs. 
5-8. 

RecentljEI we have proposed a unitarity condition method for constraining the 
CKM matrix entries, method that fully exploits the unitarity properties of CKM 
matrix, it is an exact one, its use does not depend on any approximation, and can 
be applied to both quark and fermion sectors. 
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The main goal of the paper is to make a comparison between the two approaches 
alluded to the title, and for doing that we reformulate the standard unitarity triangle 
method by requiring that the triangles sides have to depend on four measurable 
quantities. The measurability property was first raised by Jarlskog 10-12 who has 
shown that in the 3-flavors case only the functions that depend on the moduli | Uij | 
and/or the "angular looking objects" U a jUf3kU* k Upj are re-phasing invariant and 
measurable. We remark that the common four independent parameters used in the 
standard unitarity triangle approach are |?7 US |, |t/ u b|, \U c b\, and the phase 8, see 
Refs. 6-8 and the references cited therein, where the notations are the well-known 
ones. Concerning S it is not seen as being directly measurable, at least in a single- 
type experiment, even it en ters in a measurable quantity, the Jarlskog invariant 
J = Im(U a jUf3kU* k Upj^^. The J expression suggests that the "natural" four 
independent parameters to be used in any global fit have to be the moduli. Thus 
in the following we will always use as parameters four independent moduli \Uij\, 
and in particular we will provide another measurable expression for 6. That choice 
allows us to find a simple criterion for the separation of double stochastic matrices 
from the unistochastic ones, problem that to my knowledge was not yet considered 
in the physical literature, and its fulfillment is the key point in any global fit of the 
CKM entries. 

The paper is organized as follows. In the next section we find the gauge group 
of unitary matrices, i.e. the group of elementary transformations whose action on 
the unitary matrices does not change their unitary prope rties and physical content. 
In Sec. 3, we expose the unitarity condition method I I and provide the necessary 
and sufficient conditions for discrimination between the doubly stochastic matrices 
and the unitary matrices. We show that the unitarity condition is a very strong one 
and doing a global fit is not an easy matter. In Sec. 4, we reformulate the standard 
unitarity triangle approach 5-8 by using in all formulas only measurable quantities. 
In Sec. 5, we present a comparison between the two approaches, and we conclude 
in Sec. 6. 

2. Unitary matrices and their gauge group 

The n-generations CKM matrix is assumed to be unitary, and it depends on n 2 
parameters which are usually taken as n(n— 1)/2 angles and n(n+X)/2 phases, each 
set taking values in [0,7r/2], and respectively, [0, 27r). Very soorM^ it was realized 
that by a redefinition of quark fields, the elements of the first row 3x3 unitary matrix 
can be made positive. That observation was formalized as follows: the "encoded" 
physics in the CKM unitary matrix is invariant to multiplication at right and/or 
left by diagonal phase matrices D ^diagona^e 1 ^ 1 , . . . , e 1 ^'*), see e.g. Refs. 10, 15 
and 16, because we can redefine the quark fields phases. That means that we have 
at our disposal In — 1 phases whose values can be chosen by us, and the common 
choice is and/or 7r. For a non-standard choice in the 3-dimensional case see Ref 
17. That leads to n(ii — l)/2 — (2n — 1) = {n — l)(n — 2)/2 independent phases 
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and n(n — l)/2 independent angles that parametrize a unitary matrix. Here we rise 
a novel problem: are there other elementary transformations of unitary matrices 
under which their physical properties do not change?. 

Besides the left and/or right multiplication by diagonal phase matrices there is 
another transformation: multiplication at left and/or right by permutation matrices. 
Permutation matrices are matrices whose elements on each row and each column are 
zero, but one that equals unity. They interchange columns and, respectively, rows 
between themselves. Both the diagonal phase matrices and permutation matrices are 
subgroups of unitary matrices. If D is a diagonal phase matrix and P a permutation 
matrix then 

DD* = PP* = I n 

where * denotes the adjoint, and /„ is the n-dimensional unit matrix. Another 
equivalent unitary matrix can be obtained by taking the transpose of the original 
one, and if we would require that both the matrices should have similar forms 
after such a transformation, the Kobayashi-Maskawa forrrll^l will be the compulsory 
choice, and not the usual one^( We can also apply the complex conjugation to all its 
entries, and all its physical properties do not change. If we denote by T the transpose 
and by C the complex conjugation, these transformations form a subgroup of each 
other because 

T 2 = C 2 = Identity 

In conclusion the product group 

G = DxPxTxC (1) 

is the gauge invariance group of unitary matrices, i.e. the physical content remains 
invariant by applying a number of the above transformations to a given unitary 
matrix. The C-invariance has an important consequence: the range of all the (n — 
l)(n — 2)/2 independent phases is [0,7r]. 



3. Unitarity condition method 

In the following we make use of the rephasing invariance property and multiply 
the recommended form for CKM matrix™^ at left and, respectively, at right by the 
diagonal phase matrices Di — (1, e~ l& , e~ lS ) and D r = (1, 1, e %s ), and obtain the 
form that will used in the following 

(C12C13 C13S12 S13 \ 

-c 2 3S 12 e~ lS - C12S23S13 c 12 C23e~ zS - S12S23S13 S23C13 (2) 
Sl2S23e~ l<5 - C12C23S13 -C 12 S23e~ lS - S12C23S13 C23C13 / 

with Cij = cos 9ij and sy = sinSy for the generation labels ij — 12, 13, 23, and S is 
the phase that encodes the breaking of the CP-invariance. The above form depends 
on four independent phases, the phases of U21, U22, U31, C/32, instead of fives in the 
usual form. These four phases can be considered as four fundamental frequencies, 
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and their knowledge is equivalent to the knowledge of all the U entries, see e.g. Ref. 
21. 

The above model has to be supplemented by experimental data that are usually 
supplied by experimenters under the form of a positive entries matrix 



(V 2 A V 2 V 2 

v iid v us v i 



V = 



V? d V 2 V 2 b I (3) 
V t 2 d V 2 V 2 , 



More generally, the experimental data are given in terms of some functions 
fk(Vij), k — 1,...,N, that depend on the V entries. For defining a theoretical 
model we need also a linking relation between the above two objects, and we pro- 
pose 

V=\U\ 2 (4) 
The previous equation is equivalent to the following relations 



v ud 


- r 2 r 2 V 2 


_ 2 2 y2 _ 2 
— i 12 L 13' v ub — *13 




- s 2 r 2 V 2 

— *23 c 13' v tb 


2 2 
— C 13 C 23' 




v c 2 d 


2 2,2 
— s 12 c 23 ' S 13 


2 2, 
S 23 C 12 + 


2S12S13S23C12C23 cos 5, 


V 2 

cs 


2 2,2 
— c 12 c 23 "+" s 12 


2 2 
S 13 S 23 


2S12S13S23C12C23 cos 5, 


V 2 


2 2 2 , 
— S 13 C 12 C 23 + 


2 2 
s 12 s 23 


2S12S13S23C12C23 cos 5, 




2 2 2, 
— S 12 S 13 C 2 3 + 


2 2, 
C 12 S 23 + 


2S12S13S23C12C23 cos 5 



(5) 

The main theoretical problem is to see if from a matrix as JHJl one can reconstruct 
a unitary matrix as (J2J) . For that we make use of the unitarity relation 

UU* = U*U = I 3 (6) 

to obtain from the relations (JSJ the weakest form of unitarity 

J2 V 2 -1 = 0, j=u,c,t 

i=d,s,b 

E ^- 1 = > J = d,s,b (7) 

i— u,c.t 

For the moment we assume that the entries of V are such that the relations 
are exactly satisfied; then the set (JJJ is known as double stochastic matrices, and 
the subset V 2 - — \Uij\ 2 is known as unistochastic matrices**!. The constraints Q 
are necessary but not sufficient for unitarity, where from we infer that we need a 
separation criterion between the two sets. 

Furthermore, for simplicity we will choose four independent parameters to work 
with, and our choice is 

\U US \ = a, \U ub \ = b, \U cd \ = d, \U cb \ = c (8) 



I 
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Using the exact doubly stochasticity of V we form the real matrix 



S 



/ VI - a 2 - b 2 a b \ 

d v 7 1 — c 2 — d 2 c 



(9) 



■ \/a 2 +b 2 - d 2 \Jc 2 + d 2 - a 2 \J\-b 2 - c 2 j 



and observe that \S\ 2 is double stochastic. 

Making the identification Sij = Vij in relation JSJ we find expressions for Sy 
and cos *^f, the last one being given by 

r d 2 (l - 6 2 ) 2 - a 2 (l - b 2 - c 2 ) - b 2 c 2 (l -a 2 - b 2 ) 

coso = , (10) 

2a&cVl -a 2 - &V1 - b 2 - c 2 

showing that cos<5 is a measurable quantity being expressed in terms of four moduli. 
The separation criterion between double stochastic and unistochastic matrices is 
given by the physical condition 

-1<cos£<1, or cos 2 5 <1 (11) 

that is equivalent to 

a 4 - 2a 4 c 2 - 2a 2 b 2 c 2 + a 4 c 4 + 2a 2 b 2 c 4 + b 4 c 4 - 2a 2 d 2 + 2a 2 b 2 d 2 
+2a 2 c 2 d 2 - 2b 2 c 2 d 2 + 2a 2 b 2 c 2 d 2 + 2b 4 c 2 d 2 + d 4 - 2b 2 d 4 + b 4 d 4 < { ' 

relation that describes the physically admissible region in the 4-dimensional space 
generated by the moduli a, b, c, d. This result shows that four independent moduli do 
not always determine a unitary matrix, even when the moduli satisfy the relations 
0; this happens then and only then, when the relation Ulljl is satisfied. The above 
relation separates the unitary matrices from the doubly stochastic ones in our choice 
of independent moduli. The form l|10|) changes when we change the choice (jHJl. 
Relations that have to be satisfied for fulfilling unitarity have also been found in 
Ref. 23 but unfortunately they were not taken into account by the CKM fitter 
community. 

Relations Q being exactly satisfied we can compute cost) and check if the 
physical condition is satisfied. If it is, then the experimental data are compatible 
with the form (J2J and we can easily reconstruct it from the data (j3J). If the physical 
condition is violated there is no compatibility and the story ends here. 

However, the experimental data are not known with infinite precision to see if 
the relations J3J are exactly satisfied. The data being affected by errors, what could 
we do in that situation? There is the place where the gauge invariance group of 
unitary matrices enters the game. The gauge invariance group tell us that there 
is a "democracy" within the CKM matrix entries, all of them enjoying the same 
"rights" , and by consequence our choice freedom does not amount to the choice 
It is a combinatorial problem to find all the four independent moduli groups, 
and in fact there are only 58 such groups that lead to 165 different expressions for 
cos 5. Depending on the explicit group one gets one, two, three or four different 
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expressions for cos<5. The huge number for cos<5 is a consequence of the data errors 
and comes from the fact that neither relations Q are exactly satisfied, nor cos 5 
values are equal, nor in the physical region; all these three conditions have to be 
simultaneously satisfied in order to find a unitary matrix. If the data are exactly 
known and satisfy the relations JJJ) all the 165 expressions give a single number. For 
example with the choice of rational numbers 

1 1 23 26 28 

a = — = 0.22, b = = 0.00368, c = = 0.0416, d = = 0.224 

50 6250 625 125 

that are almost identical with the values recommended i J2BI the S 2 matrix is 
doubly stochastic and one gets 

7748742441690187 

cos 5 = , w 28.32 

7187500V144947356493366 

Although two conditions, the double stochasticity property and the equality of cos 6 
are exactly satisfied, the physical condition l|ll|) is not, so the above numbers are 
far away from numbers coming from a unitary matrix. If we change only one value 
and take a — 28/125 = 0.224, we find cos<5 w 1.2897. This simple test could say 
that the recommended value a — 0.22 is highly improbable to be true and has to 
be changed. However we have to wait and see what the experiments will tell us! 

To see that the condition (| 1 1 ft is constraining, we did a Monte Carlo simulation 
for Jul using as input the PDG data a = 0.2200 ± 0.0026, b = (3.67 ± 0.47) x 
10~ 3 , c = (41.3 ± 1.5) x 10~ 3 , d = 0.224 ± 0.012, i.e. the values recommended 
in Ref. 20, and the results are shown in Fig.l. The results show that only a tiny 
fraction of simulated data are in the physical region and, by consequence, a global 
fit is not easy to be done. For similar numerical results see Ref. 24. 

In the case of approximate data the explicit form for cosi5 depends on the 
independent four parameters we choose to parameterize the data and also on 
the parameters entering the last four equations in JSJ that provide in general 
different expressions for cos S, i.e. the numerical values are in general different 
cos 8^ 7^ cos 5^ , i 7^ j. In that situation we have to provide a x 2— test that has to 
take into account the double stochasticity property expressed by the conditions @ 
and the fact that in general the data are such that cos S may take different numerical 
values even for those cos i5 coming from a definite choice of the four independent 
parameters, and more, the numerical values could be outside the physical range 
(—1,1). My proposal is 



X?=£(cos^>-cos^)) 2 + £ £ ^-1 

i<j j=u,c,t \i=d,s,b 

+ E f E ' -1<co S *«<1 (13) 

j—d,s,b \i=u,c,t J 



I 



I 



2, 2008 3:24 WSPC/INSTRUCTION FILE mpla'2 



Unitarity condition method versus standard unitarity triangle approach 7 




-0.5 0.5 1 -400 -200 200 400 

cos6 cos6 



Fig. 1. Number of Events as function of cos<5. The input values are: a = 0.2200 ± 0.0026, 6 = 
(3.67 ± 0.47) X I0~ 3 , c = (41.3 ± 1.5) X 10~ 3 , d = 0.224 ± 0.012. (a) On the left panel one sees 
the number of events for the physical cos S £ (—1, 1). On the right panel are all the events and 
the interval of variation is cos5 £ [—350, 450]. The number of physical events is 1% of the total 
simulated data, and there is a peak around cos S m 0. 

We stress that the relation \\ = holds true for both double stochastic and unitary 
matrices, and only the physical condition — 1 < cos£W < i discriminates between 
them such that it has to be compulsory satisfied. The second component of the 
X 2 -test has the form 

xl = ^^_Aj (14) 

where di are theoretical functions depending on Vki, di are the measured exper- 
imental data for di, and a is the vector of errors associated to di. The formula 
X 2 = Xi + xi i s our phenomenological tool for analyzing the experimental data. 

4. Standard unitarity triangle approach 

Unitarity triangle approach exploits another consequence of the unitarity property 
Eq.©, namely the six orthogonality relations of rows and, respectively, columns 
of a unitary matrix 6 ^ 8 . Its root is in the pioneering work by Wolfensteir0 where 
he proposed a relationship between the theoretical object U and the experimental 
data available then, its nice feature being the possibility to estimate the order of 
magnitude for all U entries at a time when the experimental data were very scarce. 
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Usually one considers only the orthogonality of the first and the third columns of 
U , relation which is written as 

u ud u* ub + u cd u* cb + u td u; b = o (is) 

The main reason for that choice is the structure that arise in (|15fl by applying 
the Wolfenstein parametrization and keeping only the leading non- vanishing terms. 
However, there is another orthogonality relation that has the same structure, that 
given by the orthogonality of first and third rows 

u ud u* td + u us u; s + u ub u; b = o (16) 

that was never used in phcnomcnological analyzes. The above relations can be 
represented as triangles in the complex plane. The other four relations are neglected 
because the triangles have a flattened form in the complex plane. To see how the 
relation Ijl5(l has been exploited till now, see Refs. 7, 8, 25 and 26, and the references 
cited therein. 

The relations (|15JI - H16|I are scaled by dividing them through the middle term 
such that the length of one side is 1. The other sides have the lengths 



Ru = 



Rt 



u udUl b 

UcdU* cb 
u td u t 



bVl -a 2 -b 2 



tb 



u cd u* 



cd 

-b 2 -d 2 Vl-b 2 - 



c d 



(17) 



for Ea. lfT3|l . and respectively 



R„ 



-hi 



R, 



U u bUt b 



u us u* 



VI - a 2 - b 2 Va 2 + b 2 -d 2 

a v? 
bVl -b 2 - 



d 2 -a 2 



- Vc 2 



d 2 - a 2 ^ 

for Ea. l|16|l . On the right hand side we have written the R- values in our choice for 
independent parameters (JSJ) by supposing that the relations hold. In this case 
we can compute the numerical values on the right hand side and see if they form 
a triangle. If they do, one computes their angles and the above triangle provide us 
two independent re-phasing invariant angles. If the computed values do not form a 
triangle, i.e. they do not satisfy the inequality 



\Ru-Rt\<i<Rv 



Rt 



(19) 



the data are not compatible to the existence of a unitary matrix. This inequality 
is the analogous condition to the physical condition (fill) . Monte Carlo simulations 
for Ru and R 2 are shown in Fig. 2. They tell us that there are no problems with 
R u € [0.2, 0.6]. Concerning R 2 , it takes values in the range [—250, 250] and only 
37.5% are positive; trying to satisfy the inequalities (|19l) one finds that only 1.5% 
of them are in the physical region. The difficulties come from y/a 2 + b 2 — d 2 and 
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Vc 2 + d 2 — a 2 , that lead easily to imaginary values because a ~ d, and from the 
smallness of a, b, c and d. In the usual approach they are approximated by \Utd\, 
and, respectively, by \U t8 \, because a 2 + b 2 ~ d 2 = 1 — \ U ul j\ 2 — d 2 ~ \U t d\ 2 , etc., such 
that the imaginary values are never seen. 



2250 - 




4 00 



Fig. 2. Number of Events as function of R u , and R^. The input values are: a = 0.2200 ± 
0.0026, b = (3.67 ± 0.47) X 10" 3 , c = (41.3 ± 1.5) X 10" 3 , d = 0.224 ± 0.012. R u takes values 
within the interval R u e [0,2, 0.6], and R% e [-250, 250], and only 37.5% of them are positive. 
The condition 1191 reduces the number of Rt events to about 1.5%. 



Similar to the unitarity condition method, when the data are not exactly known, 
we have to calculate all the Ri forms corresponding to the 58 groups of four inde- 
pendent moduli. With the choice 

\U ud \ =e, \U us \=a, \U cd \=d, \U cs \=f 

one gets 



n-u — 



d^l-d 2 -f 2 



R = VI - g - rfV« 2 + d 2 + e 2 + / 2 ~T 
d^l-d 2 -f 2 

which is quite different from H17|) . so we have to require R$ « Ru\ i ^ j and 
respectively r[ 1 ' « Rf\ i 7^ j- Similar to the preceding case the problem is caused 
by \/\ — d 2 — f 2 that takes imaginary values because / is quite large, and the 
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determination of the two ratios 1(20(1 is difficult. In this case of standard unitarity 
triangle one can define a % 2 — test as follows 

x! = E(^ ) -^ ) ) 2 + (4 i) -i4 i) ) 2 + E ( E ^ - 1 

i<j j—u,c,t \i=d,s,b 

, \R$-B? ) \<1<rU+^ (21) 




Also we have to consider all the six orthogonality relations on the same footing 
and we have to include this information in the form J2J. The necessity for doing 
that comes from the fact that one orthogonality relation provides only two rephasing 
invariant angles. Or it is well knowrl^that there are necessary four such angles for 
a complete determination of a unitary matrix. Unfortunately the beautiful papeJ^ 
had no echo among the people working on global fits to the CKM matrix entries. 

Similar to the preceding method X 2 = xi + xl can De use< i as a X 2 ~test for the 
standard unitarity triangles approach. In principle both the approaches, unitarity 
condition method and unitarity triangles method, have to give the same results 
modulo small errors. Since some fitting methods could provide numerical values for 
all Vij, the unitarity test — 1 < cos 5 < 1 has to be checked. 

5. Comparison of the two methods 

Doing a global fit for the CKM matrix entries is not easy at all. The difficulty comes 
from the formulas 1)10(1 . I(17|l . ((18(l . ((20(1 and those similar to, that lead to unphys- 
ical values and we have to find numerically safe methods to satisfy the physical 
constraints 1(11(1 and/or 1(1 9(1 . Concerning S the Monte Carlo simulations, see Fig.l, 
show that the overwhelming majority of its values lies outside the physical region 
((11(1 . For other choices of the independent parameters, cos 6 can even take imagi- 
nary values as happens when on computes it by using central \Uij \ values published 
m some fitP. A similar situation arises for R u and/or R t . Fig. 2 shows that i? 2 
takes negative values in 62.5% cases, i.e. Rt gets imaginary. That happens because 
a ss d. If we try to fulfill also the condition ((19(1 then the number of physical events 
is around 1.5%, i.e. of the same order of magnitude as for cos<5. Perhaps it is a 
matter of taste to consider that it is numerically easier to satisfy the relation ((11(1 
and more difficult to implement the inequalities 1(19(1 : perhaps it comes from my 
dissatisfaction concerning the current use of the unitarity triangles approach. 

The main reason for looking for an alternative approach to global CKM fit- 
ting techniques is that there is no simple relation between the parameters defined 
by Wolfensteirpl A, A, p, 77, and directly measurable quantities. The weakness of 
Wolfenstein parametrization was intuited many years ago by Branco and LavourE^I, 
but, unfortunately, it was not fully exploited, although they gave a Wolfenstein-type 
parametrization with four moduli!. Their proposal was to choose another parame- 
ter, q = \U t d/{U us U c b)\ 2 , instead of the usual one 5 — arg([/ n f, U cs U* s U* b ), the first 
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being more relevant to B— mesons physics and a measurable quantity ^ . Of course, 
in principle, S is measurable, see the relation ()10(l that expresses cos S through mea- 
surable quantities, but until now it cannot be measured in a single experiment; we 
have to use data from many different experiments t o ob tain its value. A similar 
situation is that of measuring the Jarlskog invariant J^that depen ds on s in 6. 

As concerns the nowadays use of the unitarity triangle methocP^EEl i n mv 
opinion, its main drawback are the approximations V a 2 + b 2 — d 2 w \Utd\ an d 



d 2 



\Uts\, and similar to, for other choices of the independent parame- 



ters. In some happy situation as for example when the range of one R u is positive, 
see e.g. the case <|17|) . one can "control" Rt to obtain reasonable values for some 
physical quantities, even if R t can take imaginary values. However the problem is 
not simple and obtaining a physical result requires a very precise matching between 
the involved parameters. If we look at the triangle inequality l|19[) we easily see that 
the equality a = d is excluded since R u + Rt < 1, for a, 6, and c around the measured 
values. A physical solution is possible if their difference a — d ~ 10~ 4 — 10 -5 , i.e. a 
and d are strongly correlated taking values in a small window; se for a similar 
result. On the other hand the global CKM fits practitioners, see e.g. Refs. 28-30, 
use expressions as 

\2 



Ru 



Rt 



UudU* b 


= (i 


A 


UcdU* cb 




u td u; b 


i 


v td 


UcdU* b 


~ A 


Vcb 



Vcb 



(22) 
(23) 



that are not sensitive to such small variations of a and d parameters. In the second 
case (|18|) both R can take imaginary values and they cannot be "guessed" with 
enough precision for a reasonable fit, and we think that this is the real reason 
why the second unitarity triangle is not used, although both of them have the 
same magnitude in / i l l However, since one unitarity triangle provides only two 
invariant angles, we need to use at least two such triangles for the determinations 
of oth er t wo invariant angles for a complete determination of the CKM unitary 
matriJ^l. 



6. Conclusions 

In the paper we have shown that the unitarity condition methocP^^and the unitar- 
ity triangles method 6-8 are both a consequence of the assumed unitarity property 
for the CKM matrix. In the same time we have reformulated the second method 
such that the parameters Ri should depend on four measurable quantities, namely 
four independent |?7,y|. 

We have shown that in the real case of approximate experimental data we have 
to use the gauge invariance property of unitary matrices for imposing all possible 
constraints implied by the unitarity property. For example the parameterizations 
U13 = sise~ lS and U13 — si2S23e~ tS — C12C23S13 have to be treated on the same 
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footing, and have to lead to the same physics. 

We provided % 2— tests for both the methods under the form (|13|1 and 121JI . 
forms that are imposed by the errors that affect all the experimental data, and by 
the gauge invariance group of unitary matrices. 

By Monte Carlo simulations we have shown that the unitarity property requires 
a very fine tuning between all the CKM matrix elements, and our proposed test for 
unitarity implies the fulfillment of two conditions \i ~ and — 1 < cos S < 1. 

In our formulation both the methods are exact, they do not depend on any 
assumptions concerning the smallness of some parameters and can be used for global 
fits in both the quark and lepton sectors. 

Last but not least we made a clear distinction between the theoretical model 
given by the matrix U, and the experimental data matrix V , problem that is gen- 
erally overlooked. The unitarity is assumed and built in the parametrization of the 
CKM matrix U, but the results of the fits could violate the unitarity as the Monte 
Carlo simulations and the numerical examples in section 3 prove; see also^l where 
the central values published i rP led to a few imaginary values for cos d. In this sense 
the most important result for global fits of the CKM matrix entries is the separation 
criterion between the double stochastic matrices and unistochastic ones, problem 
that, to my knowledge, was not considered until now in the physical literature, so 
much the less by the CKM fitting community. 

Acknowledgments 

The author would like to thank C. Alexa and S. Constantinescu for a quick and 
pedestrian introduction to Monte Carlo simulations, we also thank K. K. Phua for 
inviting me to submit this contribution to Modern Physics Letters A. 

References 

1. M. Battaglia, A. J. Buras, P. Gambino, and A. Stocchi (Eds.), Proceedings of the 
Workshop The CKM Matrix an d the Unitarity Triangle, 13-16 February (2002), 
CERN, Geneva, hep-ph/0304132 

2. H. Abele and D. Mund (Eds.), Proceedings of the Two-Day-Workshop Quark-Mixing, 
CKM- Unitarity, September 19-20 (2002), Heidelberg, hep ph/0312121 

3. P. Ball, J. M. Flynn, P. Kluit, and A. Stocchi (Eds.), Proceedings of the 2nd Workshop 
on the CKM Unitarit y Triangle, IPPP Durham, A pril 2003 

4. CKM2005 Workshop, fhttp: / /ckm2005.ucsd.edu"71 

5. L. Wolfenstein, Phys.Rev.Lett. 51 1945 (1983) 

6. A. J. Buras, M. E. Lautenbacher and G. Ostermaier, Phys.Rev D50 (1994) 3433 

7. M. Ciuchini et al., JHEP 0701 013 (2001) 

8. A. Hoc ker, H. Lacker, S. Laplace, and F. R. Le Diberder, Eur.Phys.J. C21 225 (2001) 

9. P Dn,a, hep-ph/0408013 

10. C. Jarlskog, Z. Phys. C 9 29 491 (1985); Phys.Rev.Lett. 55 1039 (1985) 

11. C. Jarlskog, Phys.Rev. D 35 1685 (1987) 

12. J. D. Bjorken and I Dunietz, Phys.Rev. D 36 2109 (1987) 

13. P. Dita, hep-ph/0502125 



2, 2008 3:24 WSPC/INSTRUCTION FILE mpla'2 



Unitarity condition method versus standard unitarity triangle approach 13 

14. L. Maiani, in Proceedings of the International Symposium on Lepton and Photon In- 
teractions at High Energies, Hamburg 1977, DESY, Hamburg, 1977 

15. I. Dunietz, O. W. Greenberg, and D. Wu, Phys.Rev.Lett. 55 2935 (1985) 

16. J. F. Nieves and P. B. Pal, Phys.Rev. D 36 315 (1987) 

17. F. J. Bottela, G. C. Branco, M. Nebot, and M. N. Rebelo, Nucl.Phys. B 651 174 
(2003) 

18. M. Kobayashi and T. Maskawa, Progr.Theor.Phys. 49 652 (1973) 

19. L. L. Chau and W. Y. Keung, Phys.Rev.Lett. 53 1802 (1984) 

20. S. Eidelman et al, Phys.Lett. B 592 1 (2004) 

21. R. Aleksan, B. Kayser and D. London, Phys.Rev.Lett. 73 18 (1994) 

22. A. W. Marshall and I. Olkin, Inequalities: Theory of Majorization and Its Applications, 
(Academic Press, New York, 1979), Chapter 2 

23. G. C. Branco and L. Lavoura, Phys .Lett. B 208 123 (1988) 

24. F. J. Botella et al, |hep-ph/0502133| 

25. J. Charles et al, (The CKM Fitter Group), hep-ph/0406184 

26. M. Bona et al, UT fit Collaboration, hep-ph/0501199 
http: /www. utfit.org 

27. G. C. Branco and L. LavouraRhysRev. D 38 2295 (1988) 

28. A. J. Bura s, F. Parodi and A . Stocchi, JHEP 0301 (2003) 029 

29. A. Stocchi, [hep^ph70405038 

30. R. Fleischer, hep-ph/0405091 



